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Abstract 

We consider different configurations of ac driven quantum dots coupled to superconductor leads where Majorana fermions can exist 
as collective quasiparticles. The main goal is to tune the existence, localization and properties of these zero energy quasiparticles 
by means of periodically driven external gates. In particular, we analyze the relevance of the system and driving symmetry. We 
predict the existence of different sweet spots with Floquet Majorana fermions in configurations where they are not present in the 
undriven system. 
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1. Introduction 

There are condensed matter systems which can hold collec¬ 
tive quasiparticles that are their own antiparticles, therefore sat¬ 
isfying the Majorana condition IHIEIEI. These quasiparticles 
are termed Majorana Fermions (MFs) and follow non-abelian 
statistics. Detection of MFs in solid state systems have been re¬ 
cently experimentally proposed ElISlEl. Recently, the interest 
in encoding a qubit in these kind of excitations has grown due 
to the possibility to be non-local, a property which has a great 
potential in quantum computation due to the robustness of the 
qubit against local perturbations Q. Furthermore, how to tune 
MFs in condensed matter systems is one of the main purposes 
of research in the emergent field of topological quantum com¬ 
putation. 

In the last years, different works have shown how the appli¬ 
cation of ac fields enriches the properties of these quasiparti¬ 
cles and facilitate their tunability. For instance, it is possible 
to generate Floquet Majorana fermions (FMFs) as steady-states 
of non-equilibrium systems which present interesting properties 
for quantum computation: non-locality and non-abelian statis¬ 
tics |[8l[9l. 

In every system with particle-hole symmetry, the quasiparti¬ 
cles come in pairs = yp, therefore they can hold MFs as 
long as the energy can be tuned to zero. One of the simplest 
and most tunable system with particle hole symmetry is a dou¬ 
ble quantum dot (QD) connected via an s-wave superconductor 
Col. It is well known that the proximity effect induces Cooper 
pairs correlations across the dots CIIIEI generating effectively 
superconductivity HU. Interestingly, Fractional Josephson ef¬ 
fect, a signature of the presence of MFs |[5]|T4l|T5l, in a quadru¬ 
ple quantum dot in the presence of an s-wave superconductor 
has been predicted by Markus Biittiker and coworkers ca. 

The advantage that configurations of a few QDs connected to 
s-wave superconductors present, in order to generate and detect 
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Figure 1: Scheme of three QDs coupled by tunnel and coupled to a super¬ 
conductor. The existence of Cooper pairs generates correlations of the type 
di,o-di+i,d- in the effective Hamiltonian for the QDs. The applied magnetic fields 
and their directions are also shown in the picture. The angle 6 controls the ratio 
(see text below). 


MFs, in comparison with nano-wires ifTTl [TSl [T^ or long QD 
chains ||20| |22l proposals is their great tunability, while in 

the latter the MFs have topological protection. 

In this paper we analyze two different configurations of QDs 
in proximity to superconducting leads such that Cooper pair 
correlations are induced between the neighboring dots as long 
as the coherence length is larger than the distance between 
them. We include periodically driven gates and search for the 
conditions for appearance of FMFs. The paper is organized as 
follows: In section]^ we present the model, in sectionwe 
discuss the generation of FMFs in a double and a triple super¬ 
conducting QD. Finally, we present our conclusions in section 

IS 

2. Undriven system 

Systems of QDs coupled to s-wave superconductors have 
been a subject of study iniiniiisi because the proximity effect 
induces Cooper pair correlations that can be easily detected due 
to the low number of degrees of freedom in QDs. In a system 
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where neighboring QDs are coupled through superconducting 
reservoirs as in Fig[^ in the limit of large superconducting gap 
the superconductors can be traced out and an effective Hamil¬ 
tonian for the dots is obtained 12311241 

H = ( 1 ) 

i,(T 

^ "I" h.C.^ , 

i,cr 

which already contains effective superconductivity between 
neighboring dots. The fermionic operator represents the 
annihilation of an electron in the i-QD with spin cr. The sym¬ 
bol d- means the opposite spin to cr, which can be cr /li 

is the onsite energy in i-QD, the parameter tij+i is the effective 
tunneling probability from dot i to dot i 1 through the super¬ 
conductor by virtual occupation of the above gap excitations 
and A//+! is the effective superconducting amplitude due to the 
superconductor connecting the i and i -h 1 dots. If a large mag¬ 
netic field is applied to the dots only one spin comes into play. 
However, the magnetic fields have to be non-collinear in order 
to have s-wave type Cooper pair correlations (see Fig(TJ OOl. 
In this configuration, it is more natural to work in the basis of 
the quantization axes given by the magnetic field in each dot. 
For that purpose, we have to perform the rotation 

6 6 

d 2 ,o- cos -d 2 ,o- + cr sin -d 2 ,^ (2) 

as the magnetic field in the central QD forms an angle 6 with 
the magnetic fields in the left and right QDs (see Fig[^. The 
low-frequency hamiltonian will be given by eq. Q by neglect¬ 
ing the contribution from the high-energy spin direction in each 
dot (keeping cr =|): 

H = YjfiidJdt (3) 

i 

+ X 4m d^ df+i A^- -l- /z.c.j , 

i 

where J,- = cos f and A'= A,y+isinf. 

Therefore the normal and superconducting tunneling ampli¬ 
tudes are renormalized and their renormalization depends on 
the angle between the magnetic field directions. This depen¬ 
dence introduces a simple way to tune externally the coupling 
parameters of the system lITOl . 

In order to obtain the excitation spectrum of the sys¬ 
tem the Hamiltonian is written in the Nambu basis 'F 
(di, d\, ..., ..., dj^, d]^ as: 

i 

For a triple QD h reads 
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The eigensystem of h (h\i = T/V/) determines the quasiparticles, 
given by yt = v/ • 'F. A zero-energy solution, 4/ = 0, implies the 
presence of a pair of Majorana quasiparticles. 

In the case of a double QD one can choose an angle such that 
A' 2 = ±t[ 2 and if yUi =0, there are two MFs given by 

n = -^{diTdl), (6) 

where 6 = Only in the case where yU 2 = 0 the MFs are 
spatially separated ca. In the case of a triple QD, assuming 
A'= 0. there are two MFs given by 

(d^ ± — Of ^d2 ± (3 {^i ± d^^ 

where a = and {3 = In the case where yU 2 or yus 

^^ 2,3 '^^ 12^23 

are zero the MFs are spatially separated (251. Interestingly, the 
manipulation of the onsite-energies allows to change the local¬ 
ization of the MFs, which would be relevant for their detection 
in transport ma. 

3. Floquet Majorana fermions 

In the following, we will apply external ac fields in order 
to change periodically the onsite energies of the QDs and in 
this way obtain FMFs as steady-state solutions of the non¬ 
equilibrium problem. 

For every system described by a time-periodic Hamiltonian a 
set of solutions exists, called Floquet states, which have the 
form \il/n{t)) = where \Unit)) are time periodic func¬ 

tions called Floquet modes and are the so called quasiener¬ 
gies |[26l|27l[28l. As the quasienergies are only defined modulo 
Q, where Q = y and T is the period of the Hamiltonian, a 
system with particle-hole symmetry (with excitations in pairs 
= je) will hold FMFs if 6 = 0, +Q/2. If the frequency is 
large enough, it is a good approximation to consider the time- 
averaged Hamiltonian to describe the dynamics. For lower fre¬ 
quencies, where multi photon processes are relevant , the dy¬ 
namics becomes more involved but there is also a way to find an 
effective time-independent hamiltonian which includes as many 
photon processes as necessary (2911^ . 

The motivation to consider periodically driven quantum sys- 
(4) terns is the fact that their time-evolution is governed by an ef¬ 
fective time-independent Hamiltonian, whose properties can be 
engineered according to the particular purposes. This method, 
called Floquet engineering, has been employed to achieve 
dynamic localization EiiEai |33l, photon-assisted tunneling 
(j^ [34l or nobel topological band structures (35] |36l [571 Ell 

^3^ EasQisii. 

The application of degenerate perturbation theory in the ex¬ 
tended Floquet Hilbert space provides a high-frequency expan¬ 
sion (in powers of for this effective Hamiltonian, such as 


ri = 

72 = 
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= 2^0 EH- With the definition of the Fourier compo¬ 
nents of the time-periodic Hamiltonian 
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( 8 ) 


the leading orders of the expansion for the effective Hamilto¬ 
nian are 
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These terms will be considered below in order to obtain FMFs 
in two different configurations of driven quantum dots: DQDs 
and TQDs. 

The time periodic perturbation applied to the i-QD is: 


V{1) = ^ A, cos (Qt + tpi) dj di. (10) 


In order to study the effect of an external driving at high fre¬ 
quency, it is convenient to move to the interaction picture which 
transfers the time-dependence to the tunneling terms by means 
of the unitary transformation: U(t) = expj-/^ V{t')dt'y Only 
the non-diagonal elements change under the transformation de¬ 
pending on whether they commute or not with the time-periodic 
term: 


\d] d,+ 1 , y(o] = d]di+ 1 (A ,+1 cos (Q/ -i- <pi+ 1 ) 

-Ai COS (Qt (fi)) , 

[didi+i, V(^)] — didi+i cos -l- (ft+i) 

-\-Ai cos (Qt (fi)) . (11) 

Therefore, the renormalization of the tunneling and the super¬ 
conducting pairing depends on the symmetry of the driving i.e., 
on the intensities applied in the different dots and on the phase 
difference of the ac gate voltages between the different dots. As 
an example of this, in the case of two QDs if an ac gate potential 
is applied to each of them with the same amplitude (Ai = A 2 ) 
and frequency, the tunneling term does not change if the phases 
are equal but it does if the phase difference is tt and the opposite 
happens for the superconducting pairing (see eq[TT]) 142114^ . 

3.1. Superconducting double QD: 

In the present work, we are interested in a configuration 
such that both the tunnel and the superconducting amplitudes 
are equaly renormalized by the ac voltages. By inspection of 
eq. O, one can see that this corresponds to driving one of the 
gates periodically, it means A\ = Aq and A 2 = 0. In this case, 
the Fourier components of the time-dependent Hamiltonian are 
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Figure 2: Quasienergy gap 6q for a superconducting double QD as a function 
of the amplitude and frequency of the driving. The dark regions corresponds to 
closed gap, i.e., zero quasienergy. The plot shows that the 4-fold degeneracy 
at high frequency at the zeros of the Bessel function Jo (^) splits into two 
different sweet spots with FMFs as the frequency decreases. The bottom plot 
shows the region around the first zero and the upper plot around the second 
zero. Parameters: = F2 = 0, A'^ = T ^12 “ the 

energies are in units of A'^ which is set to 1. 


where = 42 ^™ (^)’ *e m- 

Bessel function of first kind. The zeroth order effective Hamil¬ 
tonian only predicts spatially separated FMFs if = 0 and 
A' 2 = ±t[ 2 (see eqj^. However, the following order correc¬ 
tions allow to generate new sweet spots for FMFs. The first 
order correction is zero and the effect of the second one is the 
renormalization of t[ ^ and A' ^ to some effective values given 
by: 



where / is a function of all the Bessel functions. Consider¬ 
ing only two sidebands {J'n-n = 0 for ^ > 2) its analytical 
expression becomes: 

As the ratio between the intensity and the frequency of the ac 
field increases more terms contribute to / . 

The key point in the previous discussion is that the renormal¬ 
ization of t [2 and A '2 t>y the ac field makes it possible to 
choose the driving amplitude such that Aqq = +4ff even when 
A' 2 ^ 2 - exactly what we observe in the quasienergy 

spectrum (see Figj^. In this calculation, the on-site energies pi 
and p 2 are set to zero and the static normal and superconducting 
tunnelings are different, i.e.. A' 2 ^ 2 - frequencies all 

the quasienergies are zero at the zeros of the function 27o (^) 
(approximately ^ = 2.40,5.52,8.65...) and there are no FMFs. 
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Figure 4: Lower part of the quasienergy spectrum for a superconducting triple 
QD as a function of the amplitude of the driving. The dotted (blue) line corre¬ 
sponds to feff “ Aeff and the dashed (red) to reff + Aeff. We show that the FMFs 
appear close to the conditions Aeff = ±feff- Parameters: n\ = ni - 0, nz = 1.5, 
A = 1, ^ = 0.8, n = 4, (^3 - (^1 =0. All the energies are in units of A, which is 
set to 1. 


Figure 3: Quasienergy gap 5q for a superconducting triple QD as a function of 
the amplitude and frequency of the driving. The dark regions corresponds to 
closed gap, it means, zero quasienergy. The plot shows that the 4-fold degen¬ 
eracy at high frequency at the zero of the Bessel function Jo (^) splits into 
two different sweet spots with FMFs as the frequency decreases. The bottom 
plot shows the region around the first zero and the upper plot around the second 
zero. Parameters: m = id 2 = 0, jds = 1.5, A = I, t = 0.8, (p 3 - (fi = 0. All the 
energies are in units of A, which is set to 1. 


As the second order correction becomes important, i.e., as the 
frequency decreases, two different driving amplitudes allow for 
the condition required to the existence of FMFs: the one for 
which Aeff = feff and the one for which Agff = -fefF- This is 
why at lower frequencies there are two quasienergy gap clos¬ 
ings around each zero of the Bessel function, i.e., two different 
sweet spots (the bottom panel of Fig|^ shows the gap around 
the first zero, ~ 2.40 and the upper panel around the second one 
~ 5.52). In the following we generalize this method for gen¬ 
eration of FMFs to a largest system, i.e., to an array of three 
QDs. 

3.2. Superconducting triple QD: 

Analogously to the case of the double QD, we use the driv¬ 
ing fields such that all the non-diagonal terms of the Hamilto¬ 
nian are renormalized in the same way by the ac field at high- 
frequency. That implies driving the left and right dots with ac 
gate voltages such that Ai = A 3 = Aq and A 2 = 0. Let us choose 
for simplicity t = t'l 2 ~ ^ 2 ?, ^ = ^'\2 ~ ^ 23 - S®" 

ing to analyze the presence of FMF as a function of the different 
parameters of the present setup, in particular of the phase differ¬ 
ence between the ac voltages. With the driving fields in phase 
(fi = (fi 2 -(fix = 0, the Fourier components of the time-dependent 
Hamiltonian are 
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where tm = tj'm and A^ = A27m (^)- Due to the driving 
symmetry, if we keep only the zero order term of the expan¬ 
sion for the effective Hamiltonian all the non-diagonal terms 
vanish at the zeros of 27o so there is no effective tunneling 

or superconducting pairing and the quasienergies are ±pi for 
/ = 1,2, 3. In the following, we show how the higher order cor¬ 
rections to this high-frequency approximation generate FMFs 
around these zeros. We will focus on the case pi = p 2 = 0, 
yU 3 A 0 and A' A ±t'. such that there are not MFs in the static 
case. In Fig[^ we plot the gap of the quasienergy spectrum 
as a function of the amplitude and the frequency of the driv¬ 
ing. In the limit of high-frequency the effective tunneling and 
superconducting pairing are zero so there is a four-fold degen¬ 
eracy at 6 = 0 and there are no MFs. At lower frequencies, 
these zero-quasienergy pairs appear at different amplitudes in 
which the Majorana condition is satisfied, two different sweet 
spots. This is due to the second order correction to the effective 
Hamiltonian. The largest effect of this term is a correction of 
the tunneling amplitudes, which become: 
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One difference with the double QD system is that in this case a 
small effective tunneling between dots 1 and 3 appears due to 
virtual processes. The expression for this long-range tunneling 
is 
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(19) 


Moreover, the chemical potentials p 2 and are shifted, such 
that 
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Figure 5: Spatial location of the second FMF appearing when reff = Aeff. The 
blue, green and orange bars are respectively ai, a 2 and ^3 (see eq.{^). Param¬ 
eters: fii = jU 2 = 0, A = I, t = 0.8, Q = 4, (p 3 - (fi = 0. All the energies are in 
units of A, which is set to 1. 


The shift in the chemical potentials only changes the localiza¬ 
tion of the states (see eq|7]) and the effect of the long-range tun¬ 
neling is small. In order to probe this, we plot in Fig|^ the 
quasienergy spectrum around zero and the functions Agff ± 
as a function of the driving amplitude. The sweet spots are very 
close to the zeros of these functions, indicating that the effect of 
T 13 is small. Finally, we calculate the localization of the FMFs 
found in this configuration. We choose the FMF that appears 
when = Agff (left zero in Fig|^. The Majorana pairs are 
given by: 


* 22 ) 

72 = <^i {di - 4- a2 (d2 - 4- {d^, - , 

with normalization 2 ( 2 ^ = 1. In the bar diagram in Fig|^ 
the value of the constants Ui for different values of the chemical 
potential yUs is plotted. Interestingly, Figj^shows that for certain 
values of the chemical potential yU 3 the two FMFs are spatially 
separated and that it is possible to tune the position of 72 - 
Furthermore, as we will see below, the phase difference be¬ 
tween the local ac gate voltages within each dot plays an im¬ 
portant role. Then, in order to conclude the analysis about the 
generation of FMFs in a triple QD configuration, we will show 
that the existence of sweet spots depends on the relative phase 
between the driving fields. When the two fields have opposite 
phases, ^ 3-^1 = tt, the zero order term of the expansion (|^ 
does not change respect to the previous case where cp^-cpi = 0 . 
However, the following corrections depend on the phase differ¬ 
ence. We have calculated the effective tunneling amplitudes in 
the case ^ 3-^1 =71 and the result is 
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Therefore the functions Agff become zero for the same value 
of Aq in contrast with the previous case for ^ = 0. In Figwe 
show the gap of the quasienergies as a function of the phase 
difference and the amplitude of the driving field. The measure¬ 


^0 

n 


5.8 

I 

5.3 

4.8 

2.7 

2.2 

1.7 




—7r/2 


0 


7r/2 


<5o(10-^) 


6 

i 0 

i. 

10 

5 

0 


Figure 6: Quasienergy gap Sq for a superconducting triple QD as a function of 
the amplitude and the relative phase of the driving fields cp = - (pi. The dark 

regions correspond to closed gap, it means, zero quasienergy. The plot shows 
that one of the sweet spots varies with the phase, while the other does not. The 
bottom plot shows the region around the first zero of the function JTo (^) and 
the upper plot around the second zero. Parameters: m = H 2 - 0, nz - 1-5, 
A = 1, f = 0.8, n = 4. All the energies are in units of A, which is set to 1. 


ment of this ^-dependence would be an important signature of 
the existence of FMFs. 

The existence of these exotic dynamical quasiparticles can be 
detected by connecting two metallic leads and measuring trans¬ 
port disiiiiiiii. The signatures of FMFs will be present 
in the differential conductance measurement by the fulfillment 
of the Floquet sum rule B 6 l . It is expected that FMFs could be 
measured by transport by tuning the parameters of the ac driv¬ 
ing and therefore the normal and superconducting couplings. 


4. Conclusions 


To summarize, we have discussed the existence of FMFs in 
two different configurations of QDs driven by ac gate voltages 
and coupled through superconductor leads. The simplicity of 
these systems and their tunability in comparison with other pro¬ 
posed setups which provide MFs deserve to consider them as 
suitable solid state devices to host MFs. We have shown the 
existence of FMFs by means of the expansion of an effective 
Floquet hamiltonian in power series. By modifying the fre¬ 
quency of the driving field applied to a double QD it is possible 
to control the existence of a series of sweet spots. Moreover, we 
analyze as well a driven triple QD and we predict the existence 
of sweet spots as a function of the relative phase of the local 
drivings. This method for FMFs generation can be extended 
to chains of QDs with more than three atoms. One would ex¬ 
pect that as the number of QDs increases, the localization of the 


FMFs changes and eq.([^ would be generalized. Experimen¬ 
tally, the recent achievements in the fabrication and control of 
triple 1471 SSI SSI EOl and even quadruple semiconductor QDs 
ED, also for driven configurations 15^ . open the avenue for 
the experimental realization of hybrid configurations with su- 
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perconductor contacts where FMFs can be experimentally in¬ 
vestigated. 
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